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ABSTRACT 


In the past year we have been investigating a new algorithm 
for the derivation of column sulfur dioxide, column ozone, and 
aerosol optical depth from the TOMS HDT archived albedos . 
Simulations of the measured albedos were made using a radiative 
transfer model, and a technique developed to separate the effects 
of the aerosol optical depth from the absorption optical depth. The 
retrieval adopted obtained values for the three parameters with 
high accuracy. This retrieval was then applied to the TOMS data. 



SCIENCE ACTIVITIES 


In the past year we have been developing a new algorithm to 
obtain S0 2 from the measured TOMS radiances . The S0 2 cloud that 
results from the eruption of volcanoes is largely confined to the 
altitude range above 15 km, and thus the effect of S0 2 on the 
measured N value should be the same as if one had added extra ozone 
to the altitude profile. Our first approach to the derivation of 
S0 2 was therefore to prepare a set of N-value look-up tables using 
the TOMS equatorial ozone profiles (our test case was the El 
Chichon eruption), plus a range of column sulfur dioxide values, 
for particular solar zenith angles, scan angles, and azimuthal 
angles, etc. . An example of such a table for the A and B pairs used 
in the current algorithm is given in Table 1. 

An examination of Table 1 shows that, for each pair, there 
will be many sets of values for the ozone and sulfur dioxide column 
amounts that will lead to the same N-value. If we plot the column 
ozone amounts versus the corresponding column sulfur dioxide for a 
particular pair then the result is a straight line, at least under 
the conditions that we have studied. If we now plot the second pair 
on the same graph, the interaction of the two lines will give the 
actual column amounts for ozone and S0 2 . Figure 1 shows such a 
plot. 

One thing that is immediately clear from Figure 1 is that the 
two lines are almost parallel. The slope of the two lines is 
proportional to the ratio of the S0 2 absorption coefficient to the 
ozone absorption coefficient for the particular pair. Table 2 lists 
this ratio for the A and B pairs currently used in the TOMS 
algorithm. The ratio of the S0 2 coefficient to the ozone 
coefficient for the two pairs differs only by 3.7%, which is 
consistent with the results shown in Figure 1. We chose a modified 
set of pairs. A' and B', also shown in Table 2 which give a wider 
angle between the two lines, as shown in Figure 2, however the 
ratio of the cross sections still only differs by 15%. 

We than applied this algorithm to the TOMS HDT data for April 
15th, 1982. Figure 3 shows the results obtained for a pass over the 
volcanic cloud. Each point is the average across a scan, where the 


scan angle has been limited to less than 30 degrees. In addition, 
only data for cloud free conditions (reflectivity less than 0.15) 
has been analyzed. A significant feature of Figure 3 is that as the 
column S0 2 increases across the cloud, the column ozone decreases. 
This could be a chemical effect, but one must remember that the 
volcanic plume that brings S0 2 into the stratosphere could also 
bring ozone-poor air from the troposphere. Our concern however was 
that the effect was an artifact due to the basic assumption that by 
using pairs of wavelengths the effect of aerosols would cancel out. 

We then calculated the effect of sulfate aerosols on the 
radiances using the Dave code. Table 3 gives results for the 
retrieved column sulfur dioxide and ozone using the algorithm 
discussed above for various optical depths of aerosols . There is an 
obvious impact of the aerosols on the retrieved column amounts. We 
therefore decided to see if it were possible to account for the 
aerosols. In order to do this, we had to develop another algorithm. 

The principle behind the derivation of total ozone from the 
TOMS radiances is that most of the measured radiance comes from the 
troposphere, and not from where the bulk of the ozone resides. Thus 
all of the radiation passes through the ozone layer twice. In other 
words the efficiency for the detection of added stratosphere ozone 
is 1.0. But in essence this is the same as saying that the measured 
N value depends only on the absorption optical depth T^, in the 
stratosphere, where T^ is given by s- 

Tab, = L.(Q*a + Z*y) (!) 

where Q and 2 are the column amounts of ozone and sulfur dioxide, 
and a and y are their respective absorption coefficients . L is the 
path length. Figure 4 shows calculated N values versus the 
absorption optical depth for one particular set of solar zenith 
angle, scan angle, and azimuthal angle. In this and for all future 
calculations the profile for the tropospheric ozone has been kept 
constant. The ozone content and the sulfur dioxide content were 
varied randomly, and each point represents one of these 
combinations. As expected the points fall on the same line. Thus if 


one has a measured N value, one can determine the absorption 
optical depth. 

Sulfate aerosols are non-absorbing, and when added to the 
stratosphere only scatter radiation back to space, thus reducing 
the amount of solar flux that reaches lower altitudes, the net 
result being a decrease in the N value. The impact of the aerosols 
on the measured N value will depend on the aerosol optical 
thickness, the aerosol particle size distribution, the aerosol 
refractive index, and the altitude profile. In our initial analysis 
we have assumed that the shape of the altitude profile of the S0 2 
and the aerosols is gaussian with a peak at 25 lan, and a half width 
of 2 km. The refractive indices for the aerosols were taken from 
WMO TD #24. 

We first calculated the N value for various values of T^., and 
aerosol optical depths, T aer . Then we derived an effective 
absorption optical depth, T eff , from the pure absorption look-up 
tables similar to the curve shown in Figure 4. We then calculated 
a quantity AT, given by:- 

AT = T. ff - T^, (2) 

Figure 5 shows the calculated AT for ground reflectivities of 0 and 
0.4, a solar zenith angle of 30°, and zero scan and azimuthal 
angles for the four TOMS wavelengths. To a first approximation one 
can assume that AT is independent of wavelength, i.e. for each of 
the wavelengths one can write :- 

T.„ = T^. + A ( 3 ) 


where A is a constant. 

When we plot A versus the aerosol optical depth for a 
particular case of solar zenith angle, scan angle, and azimuthal 
angle, the curve is almost linear. We found no case for any of the 
range of angles met with in the TOMS data analysis which did not 
show a single value dependence. Having determined A from the 
measured radiances, we can now go back to calculated tables, as 



illustrated in Figure 5, to determine the exact value of AT for 
each wavelength, and hence T^. for each wavelength. It is now a 
simple job to determine the values of Q and 2. Figure 6 shows some 
preliminary results for the same conditions as Figure 3 for S0 2 , 
Ozone and aerosol optical depth. The results for S0 2 are lower than 
for Figure 3, and the corresponding dip in ozone is less. The 
aerosol optical depth curve follows the S0 2 curve. 

We then applied the new algorithm to the data for the entire 
globe. The results of this analysis are shown in Figure 7. 
Interesting features are seen over China and the Sahara, namely 
negative optical depths . Both China and the Sahara have 
tropospheric dust clouds during April and these are composed of 
silicate aerosols which are absorbing. As the assumption in the 
algorithm was that we had non-absorbing aerosols, these dust clouds 
could have an apparent negative optical depth. 



N VALUE FOR RANGE OF SO, AND OZONE COLUMN AMOUNTS 
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TABLE 3. 


INPUT OZONE : 275 D.U 
INPUT S02 : 50 D.U 


AOT=O.Q5 AOT=0.1 AOT=OJ2 


REF. SZA. SCN.AZM 
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FIGURE 5. CALCULATED AT FOR TOMS WAVELENGTHS 
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LATITUDE 

FIGURE 6 OZONE, S0 2 , AEROSOL OPTICAL DEPTH, APR. 15, 1982 
(COLUMN AMOUNT DU, OPTICAL DEPTH TIMES 100) 
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